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Abstract. A new approach to the study of the AGN phe- 
nomenon is proposed, in which the nucleus activity is re- 
lated to the metric of the inner massive black hole. The 
possibility of a Toroidal Black Hole (TBH), in contrast to 
the usual Spherical Black Hole (SBH), is discussed as a 
powerful tool in understanding AGN related phenomena, 
such as the energetics, the production of jets and the ac- 
celeration of particles, the shape of the magnetic field and 
the lifetime of nucleus activity. 
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1. Introduction 

Huge energetics of Active Galactic Nuclei (AGN) are usu- 
ally interpreted in terms of accretion around a massive 
black hole (Frank, King & Raine 1985). So far, the stan- 
dard assumption has been that the black hole is a spherical 
(SBH), static (i.e. described by the Schwarzschild metric) 
or else rotating around its axis (i.e. described by the Kerr 
metric) one. However, models of collapse leading to dif- 
ferent topologies of the resulting object have been investi- 
gated (Smith & Mann 1997) and a mathematical descrip- 
tion of their properties in the framework of general rela- 
tivity has been developed by several authors, e.g. Vanzo 
(1997), Brill, Louko & Peldan (1997). Recently, a few hints 
about the possible connection between topological models 
of black holes and high-energy astrophysical phenomena 
like AGNs have been claimed (Spivey 2000). More specif- 
ically, the possibility points toward a Toroidal Black Hole 
(TBH), in which the event horizon is topologically equiv- 
alent to a torus, to explain the AGN observations. 
In this paper, we consider such a situation and evaluate 
the influence of this topology in accounting for it. 



First, we give a basic description of the conditions needed 
to obtain a TBH and the possibilities of matching them. 
Second, the geometry of a TBH is characterized through a 
dynamical treatment of particle motion around it. Then, 
some hints about the shape and intensity of magnetic field 
are provided, yielding the emission properties of such a 
structure; along this route, a particles acceleration and 
jet production mechanism is proposed. Last, an interpre- 
tation of the finite lifetime of an AGN is deduced from the 
whole picture. 

2. The toroidal structure for a black hole. 

2.1. The spacetime metric. 

A generalization of the black hole metrics can be given as 
follows (Smith & Mann 1997): 



ds' = -iv + b-^]df 



dr 



+R^[d0^ + c sinh^{y/^e)d<j)^], 



(1) 



where i, r are the time and radial coordinates, 6 and 4> are 
coordinates on a two-surface of constant curvature and V 
is a potential term, whose meaning will be discussed later. 
The parameters 6, c, a fix the topology of the structure: in 
particular, if 5 = —a = (b = —a following the solution of 
Einstein field equations in empty space) and a — > 0, c = 
+ i, then the topology is that of a torus and the space-time 
is asymptotically anti-de Sitter (AdS). We recall that a 
Schwarzschild metric for an asymptotically flat spacetime 
exhibits 6 = and -1-1 instead of the potential term V. 
Therefore, such a configuration is bounded to the presence 
of V. A standard expression for the potential would be of 
the form: 



A , 



(2) 
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where A is a negative constant in order to provide an AdS 
spacetime. However, this cannot be the cosmological con- 



2 



F.Pompilio et al.: A toroidal black hole for the AGN phenomenon 



stant, which is known to have a positive value (Perlmut- 
tcr ct al. 1997, 1999: Ricss ct al. 1998; Rocs & Hanm- 
or-Rashid 2000) and which is a few orders of magnitude 
smaller than needed to act as an effective AdS term. That 
is why another physical reasoning must be found for V 
and it must be of astrophysical nature (e.g. gravitational 
potential of the surrounding galaxy); anyhow, some ideas 
about its origin will be further given. 

2.2. The geodesies. 

A simplified picture of the gravitational potential of a 

rotating TBH can be gained by the motion of particles 
around it. Actually, particle dynamics is the best tool 
to clarify the physical phenomena happening around the 
hole. Such a task can be pursued by solving the geodesies 

equations (Lawden 1982): 



d / dx' 
ds V ds 



dgjk dx^ dx^ 
dx^ ds ds 



0, 



(3) 



for the different components gri of the metric tensor. The 
corresponding radial equation can be stated as follows: 



dr 
di 



c\e-V) + 



+ — TT-o h^V 



1/2 



(4) 



where h is the specific angular momentum along the ^;-axis 
and fc is a constant of integration. Therefore, the energy 
conservation equation can be derived: 



2MGc^ 2MGh'^ 



which can be compared with the SBH situation: 
2MGc^ 2MGh^ 



= const. 



(5) 



(6) 



.2 

The basic difference is based on the centrifugal term ■^: 
indeed, if V has no radial dependence, then the TBH has 
no centrifugal barrier and, in order to reach the black 
hole, the accreting matter will just need to overwhelm a 
certain threshold value independent of r and fixed by the 
right-hand side of the above equation. We stress that this 
trend is of fundamental relevance for the amount of angu- 
lar momentum which must be exchanged by the infalling 
gas before accreting. 
To depict particle trajectories and in order to make a 
comparison with the SBH result, we followed the standard 
approach of parametrizing the orbit in a general elliptical 
way (Chandrasekhar 1983): 



rix) = 



(1 + e cosx) ' 



(7) 



in terms of a dimensionless parameter x along the orbit, 
then we solved the time (^) and angular (^) compo- 
nents of the geodesies, with zero initial conditions. In Fig.l 
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Fig. 1. Particle trajectories in the equatorial plane, assum- 
ing they arc experiencing just the TBH gravitational poten- 
tial; the orbits are parametrized as elliptical with p = 5.6 and 
V ~ const is Eissumed. Different values for the e parameter 
have been assumed: {solid line) e = 0.2, {short-dashed line) 
e = 1.0, {long-dashed line) e = 4.5. 



the trajectories in the equatorial plane are sketched for 
V ~ const with p = 5.6 and different choices for the e pa- 
rameter. The geodesies have been evaluated for e = 0.2, 
e = 1.0 and e = 4.5. We recall that, in general, for a 
central potential of the form V = — it holds: 



e= Wl 



2h^E 



(8) 



for a particle with energy E and mass to. Then, the first 
two choices for the e parameter correspond respectively 
to low and mildly energetic particles. In this case, parti- 
cles spiral around the torus and are focused in a circular 
region internal to the structure, whose extent is wider for 
the more energetic particles. As they turn around the torus 
they gain acceleration and, in particular, while approach- 
ing the centre of symmetry (r ^ 0), their acceleration is 
proportional to the time they spend in the inner region, 
as follows: 



dr 



s/2MG 



-3/2 



a oct 



-8/5 



(9) 



so that the less they spend in the center the more they 
can gain and we will refer to this stage as a gravitational 
kick: eventually, when a particle reaches the matching- 
condition orbit, it can be kicked out through the circular 
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inner region. We note that e = 1.0 provides a limiting 

state, because a few particles can escape the black hole 
through a parabolic orbit. Last, particles with e greater 
than unity have too much energy to get trapped by the 
black hole, therefore they just approach it and experience 
a slingshot mechanism, kicking them off. 
If they are focused to the center, particles are likely to 
undergo collisions, which could be an effective way to ex- 
change angular momentum and can provide another mech- 
anism of acceleration by means of a Fermi-type stochastic 
procedure. In that framework, the energy exchanged by a 
particle per unit time due to collisions is proportional to 
the particle velocity v and to the size I of the acceleration 
region (Longair 1981): 



dt 3 V d 



(10) 



for a particle of rest-mass mo, charge q and Lorentz factor 

7. The effect of the Lorentz force is a drift in the azimuthal 
coordinate, adding to the rounded spiral motion or even- 
tually against it and braking the flow; nevertheless, this 
motion provides a shear effect and an energy exchange 
between +/— charged flows. On the other hand, the su- 
perposition of the centrifugal force and the Coriolis force 
(opposite to the direction of rotation) can allow particles 
to move towards other orbits and could be responsible 
for particles pointing to the centre of symmetry. The az- 
imuthal drift develops a ring current associated to an en- 
circling magnetic field, which adds to the disk component. 
At the end, an axisymmetric toroidal plasma configuration 
arises, resembling the tokamak device involved in nuclear 
fusion studies. The total magnetic field will exhibit helical 
field fines (Dendy 1993): 



The conservation equation depends on the the time spent 
by the particle in that region if and if we neglect diffusion 
and source terms, it holds: 



B = W X 60 -I- B^e^,, 



(13) 



dN d , ^ N 
-^ = -J,[-ENiE)]--, 



(11) 



where the potential ijj can be inferred by means of 
the Grad-Shafranov equations, derived by the pressure- 
balance condition on the azimuthal component: 



(j X B) ■ e^p ^ Vp ■ ^ 0, 



where a= ^. The final spectrum of the particles emitted and the radial component: 
by the core in a steady state situation (^ =0) is given 
by: 



(i X B) • = — 



dp 

dr ' 



(14) 



(15) 



N{E) oc E 



(12) 



so that a jet starting from the center will have a power-law 
distribution. 

We remind that to have an effective acceleration mech- 
anism, the particles must have energies greater than the 
maximum energy loss felt during acceleration, or else the 
loss can be overcome through a sufficiently rapid initial ac- 
celeration mechanism: that is exactly what makes a TBH 
a likely candidate for an effective process. 

2.3. Extended particle dynamics. 

The spiral-shaped trajectory is not the only motion the 

particles undergo. In addition to it the centrifugal and the 
Coriolis forces must be added; moreover, if we consider 
the flow around the TBH to be embedded in the mag- 
netic field produced by the accretion disk, a Lorentz force 
term arises. Actually, the accretion disk can be considered 
to be composed of annuli of charged particles fluxes, pro- 
viding azimuthal currents: this results in a magnetic field 
in which field lines encircle those currents. Therefore, the 
three forces can be evaluated as follows: 



- Fcf = 7mo(fi xr) xil; 

- Fcor = mo 



-Fb 



2^%+r 



dt 



(Cr X n); 



{v X B); 



and the plasma is self-confined to the tokamak shape by 
means of a 6-pinch and a z-pinch effects. 
However, it's worth stressing that this tokamak-like 
plasma structure has a basic difference with respect to 
the real tokamak device, that is the strong gravitational 
field. In fact, the plasma shape is mantaincd by means of 
the TBH gravitation, in contrast to the external magnetic 
field involved in tokamaks. Therefore, a strong gravita- 
tional force is likely to prevent diffusive motion, which 
appears in magnetic confined toroidal plasma. This gives 
to the system enough stability to treat it by the Grad- 
Shafranov equations. 

In conclusion, particles spiral around the torus and eventu- 
ally approach the acceleration inner region, where they are 
kicked out by a gravitational kick or they partake in the 
Fermi-type power-law distribution and fiow away along 
the TBH axis; there, they can be trapped by the mag- 
netic field, coUimating them in a jet. 
There is another process acting in such a configuration 
related to dynamical instabilities. Indeed, particle tra- 
jectories are supposed to experience Lindblad resonances 
as follows from perturbations in the gravitational poten- 
tial (Binney & Tremaine 1987). Both with particle col- 
lisions and energy/momentum exchange, instabilities can 
act as an accretion mechanism, whenever the gravitational 
kick conditions are not matched. Gravitational instabili- 
ties could also influence the tokamak stability, but they 
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should face plasma instabilities and the competition be- 
tween the two fixes the final configuration. Actually, to 
be exhaustive and more quantitative, we should perform 
a hard magnetohydrodynamical (MHD) treatment of the 
particle motion. Since it requires a hard computational 
work, we cannot discuss it here; nevertheless, we recall 
that MHD 2D simulations have shown that turbulence in 
the accretion disc has a fundamental role. The matter in 
the region of unstable orbits does not follow simple energy 
and angular momentum conserving free fall trajectories, 
but a deep exchange is involved by means of turbulence 
on the way to the event horizon and also the shear compo- 
nents are enormous (e.g. Hawley & Krolik 2000, Igumen- 
shchcv & Abramowicz 2000 and references therein). This 
will help putting matter onto orbits which spiral around 
the torus, so that it easily reaches the inner region, where 
it will experience the gravitational kick. 

3. The emission properties of a TBH. 

The huge wide-band emission from AGN is mostly ex- 
plained through processes involving the; acc;retion disk 
(Frank, King & Rainc 1985). Nevertheless, the most en- 
ergetic emission, i.e. gamma rays, is far from being sclf- 
consistcntly explained: the observed gamma ray energy 
along the jets in quasars is close to four orders of mag- 
nitude higher than expected by theory (Spivey 2000). In 
particular, the mechanism proposed by Blandford & Zna- 
jek (1978) cannot account for such an amount of energy, so 
that it cannot be the only responsible for it. As discussed 
in the previous section, in the framework of a TBH a nat- 
ural explanation for jet production emerges, then it is im- 
portant to estimate the highest values for the Lorentz fac- 
tor that particles can reach. The limit to that value is at- 
tained because of Inverse Compton (IC) and synchrotron 
losses (Longair 1981): 



I max 

^SYNC 
I max 



\/5xl0"i^ . 
B ' 



where Urad refers to the radiation field particles move in. It 
follows that in the approximation of a keplerian orbit at a 
distance corresponding to the event horizon, the maximum 
Lorentz factor around a TBH is increased by a factor: 



..TBH 



Imax ^ jt^TBH 



Imax 



^SBH 



(16) 



for fixed values of the radiation field and the magnetic 
field. 

As concerns the spectrum, a TBH can give rise to intense 
bremsstrahlung and synchrotron emission for a sufficiently 
high angular velocity w attained by the spiralling particles. 
If we consider the aforementioned power- law distribution, 
then the emission due to relativistic bremsstrahlung, syn- 
chrotron radiation and IC will be (Rybicki & Lightman 
1979): 



- hrems OC 1/ ^ ^ ^ ; 

where B is the embedding magnetic field and vq is the 
injection energy. 

Moreover, the Blandford-Znajek mechanism is increased 
by the extra magnetic field added to the embedding disc 
field, due to the dynamics leading to the tokamak-like 
plasma configuration; as a result, the power extracted in- 
creases with respect to the SBH situation by a factor of 
(Frank, King & Raine 1985): 



Wb-z oc 



TBH 



\ BsBH J 



(17) 



In addition, a new process could be involved in emission 
properties of a TBH, due to the tokamak-like plasma, that 
is the sawteeth emission, i.e. an intense periodic burst aris- 
ing from the centre of a tokamak device, whose energy 
ranges from hard X-rays up to gamma rays (Bendy 1993). 
Its physical origin is still debated, but it seems to de- 
pend on magneto-hydrodynamical instabilities acting on 
the plasma, such that a slow increase in plasma pressure is 
followed by an abrupt fall, the current rapidly grows and 
magnetic reconnection feeds the burst. 
Another emission process, that is purely related to the 
space-time characteristics, is the Penrose process (Penrose 
1969). It follows from the fact that in a rotating black hole 
the surface on which the time component of the metric 
tensor vanishes, i.e. gu = 0, differs from the event horizon: 
the two surfaces enclose a region called ergosphere where 
the energy of a particle as observed from far way can be 
negative. The original Penrose process follows from the 
decay of a particle into two photons, one of which crosses 
the horizon and gets lost, the other one gains energy in 
the range fixed by the Wald inequality (Wald 1984): 



gttf'^ <e<^E + ^v{E^-guf'\ 



(18) 



then a TBH is more efficient than a Kerr black hole in 
producing high-energy photons in this way, provided that 

TBH > SBH 
Imax — Imax ' 

4. Lifetime of AGN activity. 

If the AGN activity is related to the toroidal shape of the 
black hole, a transition to a quiescent state is expected 
as the black hole reaches SBH status. In fact, as matter 
accretes the hole, the event horizon increases, the torus 
inflates and lately looses its starting configuration turning 
spherical. On the above line, we can infer a more quan- 
titative estimate for the lifetime of the activity phase by 
equating the metric tensor components for a TBH and a 
SBH. If Rin is an initial radial dimension for the torus. 
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e.g. the middle value of the torus thickness with respect 

to the centre of symmetry, and Rjin is the final radius 
of the SBH, then the transition happens if the following 
condition is matched: 



V- 



2MG 



or else: 



2MG 



1 



V + ^ ' 



(19) 



(20) 



where Tg is the gravitational radius. Besides, the lifetime 
is related to the accreted matter AM and to the accretion 
rate ^ by: 



dt{life) 



yielding: 



AM 

dM/dt' 



(21) 




dt{life) oc 



dM/dt 



(22) 



Since ^ is observationally estimated, we could enter the 
lifetime debate if only we knew the potential. Lacking 
that, wc can only study dt{life) as a function of some 
hypothetical functional form for it: 

1. y ~ const; 

2. V(x 1/R; 

3. V oc log R; 

where the first case could be regarded as a cosmological 
vacuum energy density, while the other expressions could 
mimic the background potential of a surrounding axisym- 
metric galaxy. In particular, the second functional form 
refers to an embedding newtonian gravitational field; on 
the contrary, the logarithmic shape is motivated in order 
to reproduce the flatness of the galactic rotation curve 
(Binney & Tremaine 1987). In this frame, the formation 
of a TBH can be deduced if a protogalaxy develops a 
sufficient extended background potential and the collapse 
to a massive black hole is forced to lead to a toroidal 
configuration^ . 

In Fig. 2, the trend of dt{life) is sketched for the above 
mentioned different potentials, assuming a constant ac- 
cretion rate dM/dt. The 1/R potential gives a fixed value 
independent of while the V ^ const situation gives 
an increasing function, as the log R potential does. If the 
AGN lifetimes have small dispersion and are strictly fo- 
cused on a single value, e.g. 10^ — 10^ yr (e.g. Cavaliere 
& Vittorini 2000), then the 1/R potential seems to give a 
likely interpretation to it. 

Last, we stress that the TBH/SBH transition is likely to 

^ We suggest that another way to obtain a seed TBH should 
follow the dynamical evolution of a black hole binaxy system; 
however, such a topic deserves a strict treatment. 



Fig. 2. AGN lifetime for a TBH model with constant accretion 
rate: {solid line) V ~ const; {short-dashed line) V oc 1/R; 
{long- dashed line) V oc log R . The normalization has been 
arbitrarily fixed and r = 



be a violent event in the galaxy evolution. From a physi- 
cal point of view, the galactic potential turns from being 
axisymmetric (TBH) and with a relevant external com- 
ponent {V) to a spherical configuration dominated by the 
central SBH. Therefore, a rapid and huge accretion should 
happen, the black hole swallows the most of the surround- 
ing, so that the activity stops and a quiescent state is 
reached. 



5. Comparison between TBH and SBH. 

In this section, we trace a few guidelines to emphasize why 
a TBH could be more convenient in explaining the AGN 
phenomenon than the standard SBH model. 
First of all, the TBH provides a powerful mean to accel- 
erate particles and to give them high Lorentz factors^; on 
the contrary, a SBH has more severe limitations in do- 
ing it and usually other complicated phenomena involving 
the surrounding galaxy and inter-stellar medium must be 
invoked, e.g. shocks induced by stellar ejections or hy- 
drodynamical expansion (e.g. Rees, Phinney, Begelman & 
Blandford 1982). Moreover, the variety of emission, which 
is well described by the standard SBH model, is not re- 
stricted by a TBH model, which affects just the energetics. 
On the other hand, our model could explain the produc- 



^ For sake of completeness, we stress that such a mechanism 
could be efficient in other astrophysical contexts, e.g. acceler- 
ation of Ultra Relativistic Cosmic Rays. 
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tion of jets in a more likely way and could account for the 
large angle dispersion observed, as follows from the open- 
ing angle of the TBH. 

An important new perspective with this model is the way 

it explains the finite lifetime of the AGN phenomenon, 
that is hardly interpreted by the standard scheme through 
a multiparameter modelling of black hole Initial Mass 
Function (IMF), accretion rate and other quantities (e.g. 
Monaco, Salucci & Danese 1999). Indeed, we stress the 
power of a TBH model in significantly simplifying the fi- 
nite lifetime puzzle. 

So far, we see no drawbacks for a TBH model with re- 
spect to a SBH model, but we underline how it could 
improve and complete the general picture of the AGN 
phenomenon. An argument which could be raised against 
our view could come from iron emission line observations, 
well explained in terms of being generated close to a spin- 
ning SBH, i.e. a Kerr black hole (e.g. Fabian et al. 1989, 
Matt, Perola & Stella 1993). Nevertheless, the emission 
line arises several gravitational radii away from the hole, 
where particles could not distinguish among a SBH and a 
TBH. However, the difference in the relativistic gravita- 
tional shift effect should be: 



and for Vg it holds: 



V - 3rg/R 



1/2 



^{TBH) oc ^^{SBH), 



(23) 



(24) 



therefore this effect is likely to be undistinguished by re- 
cent X-ray data; besides, the deep model dependence of 
the Kerr black hole interpretation of the iron line gives 
no hope for a choice among the two models, until higher 
resolution observation will be available by the next X-ray 
satellites Chandra and XMM. 

6. Conclusions. 

We inferred the main properties of a TBH in contrast to 
the widely accepted SBH model to describe the AGN phe- 
nomenon. In particular, we found that a rapidly rotating 
TBH can provide more energy to feed particle acceleration 
and production of jets. Besides, a TBH should undergo a 
transition to a SBH and so to a quiescent state (if the AGN 
activity is a pecularity of TBH alone) , forcing nuclear ac- 
tivity to fade and giving a reasonable physical explanation 
to the lifetime of AGNs. 

The only caveat of our picture is that the model relies 
on the presence of a potential V acting as an AdS term, 
whose origin cannot be due to the cosmological constant 
and is missing in the present discussion. The topic de- 
serves a more strict treatment, anyhow we tend to con- 
sider it as a term due to the surrounding gravitational 
structure, favouring the TBH collapse. Other hints could 



be found in cosmology, for instance in quintessence models 
in which the scalar field should have non minimal coupling 
to gravity, or else in cosmological models allowing dark 
matter-dark matter coupling (e.g. Perrotta, Baccigalupi 
& Matarrese 2000, Amendola 2000); nevertheless, such a 
cosmological origin is very unlikely, due to the extremely 
low value for the potential, preventing its effectiveness in 
acting as a collapse source. 
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